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Abstract: The Glycolysis Model Has Been Solved Numerically In Two Dimensions By Using Two Finite 
Differences Methods: Alternating Direction Explicit And Alternating Direction Implicit Methods (ADE And 
ADI) And We Were Found That The ADE Method Is Simpler While The ADI Method Is More Accurate. Also, 
We Found That ADE Method Is Conditionally Stable While ADI Method Is Unconditionally Stable. 
Keywords: Glycolysis Model, ADE Method, ADI Method. 



I. Introduction 

Chemical reactions are modeled by non-linear partial differential equations (PDEs) exhibiting 
travelling wave solutions. These oscillations occur due to feedback in the system either chemical feedback (such 
as autocatalysis) or temperature feedback due to a non-isothermal reaction. 

Reaction-diffusion (RD) systems arise frequently in the study of chemical and biological phenomena 
and are naturally modeled by parabolic partial differential equations (PDEs). The dynamics of RD systems has 
been the subject of intense research activity over the past decades. The reason is that RD system exhibit very 
rich dynamic behavior including periodic and quasi -periodic solutions and chaos(see, for example [8]. 

1.1. MATHEMATICAL MODEL: 

A general class of nonlinear-diffusion system is in the form 
du 

— — dihu + a x u + b x v + f(u,v) + fli(x) 
at 

dv 

d 2 Au + a 2 u + b 2 v — f(u,v) + g 2 (x) 



(1) 



dt 

with homogenous Dirchlet or Neumann boundary condition on a bounded domain £2 , n<3, with locally 
Lipschitz continuous boundary. It is well known that reaction and diffusion of chemical or biochemical species 
can produce a variety of spatial patterns. This class of reaction diffusion systems includes some significant 
pattern formation equations arising from the modeling of kinetics of chemical or biochemical reactions and from 
the biological pattern formation theory. 

In this group, the following four systems are typically important and serve as mathematical 
models in physical chemistry and in biology: 

• Brusselator model: 

% = -{b + lXh = 0,a 2 = b,b 2 = 0,f = u 2 v,g 1 = a,g 2 = 0, 
where a and b are positive constants. 

• Gray-Scott model: 

<h = -(f + k),b x = 0,a 2 = 0,b 2 = -F,f = u 2 v, gi = Q,g 2 = F, 
where F and k are positive constants. 

• Glycolysis model: 

% = -l.bt = k,a 2 = 0,b 2 = -k,f = u 2 v,g 1 = p,g 2 = 8, 

• where k, p, and 8 are positive constants 

• Schnackenberg model: 

a 1 — — k, bi — a 2 — b 2 — 0,f — u 2 v,g 1 — a,g 2 — b, 
where k, a and b are positive constants. 

Then one obtains the following system of two nonlinearly coupled reaction-diffusion equations (the Glycolysis 
model), 

d u n \ 

■ = d x Au - u + kv + u 2 v + p, (t,x) e (0, oo)xQ 



dt 
dv 

— = d 2 Av - kv - u L v + 8, (t, x) e (0, oo) x Q. 
at 

u(t,x) = v(t,x) = 0, t>0, xedQ. 
u(0,x) = u 0 {x), v(0,x) = v 0 {x), xeQ. 



(2) 



where d 1 ,d 2 ,p,k and 5 are positive constants [9]. 
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II. Materials And Methods 

2.1. Derivation Of Alternating Direction Explicit (ADE) For Glycolysis Model: 

The two dimensional Glycolysis model is given by 

du (d 2 u d 2 u\ n 

T~ = d i 3~T + )-u + Kv + u 2 v + p, (t,x) 6 (0,oo) x Q 



dt \dx 2 dy 
dv (d 2 v d 2 v , 

— = d 2 [ ~z r + ^-r - Kv - u 2 v + deZ, (t, x) £ (0, w)x£! 
at \ox A ay A 



(2) 



We consider a square region 0<x<l, 0<y<l and u, v are known at all points within and on the 
boundary of the square region. We draw lines parallel to x, y, f-axis as x=ph, y=qk, and t=nz, p,q=0, 1, 2,...,M 
and n=0, 1,2,..., N, where h=dx, k= dy, z= St. 

Denote the values of u at these mesh points by u(ph,qk,m) = u p q n and v(ph,qk,nz) = v p q n . The 

explicit finite difference representation of the Glycolysis model is 

" p,q,n+\ ~ u p,q,n d \ 2 

= ~2 < - u p+lq,n - 2u p,q,n +u p-l,q,n +u p,q+l,n ~ 2u p,q,n + u p,q-\,n > " u p,q,n + Kv p,q,n + " p,q,n v p,q,n + P 
z h 

v p,q,n+\ ~ v p,q,n d i 



Z h 

and 



2 p+l,5,n p,q,n p-l,q,n p,q+\,n p,q,n p,q-l,n> p,q,n p,q,n p,q,n 

h 



d \ z 2 
"p,q,n+l -"p,q,n = ~ p+l,q,n - 2u p,q,n + u p-l,q,n +u p,q+l,n - 2u p,q,n + " M -l,n > " zu p,q,n +zKv p,q,n + zu p,q,n v p.q.n + P 1 
h 

d 2 z 2 
v p,q,n+\ - v p,q,n = ~ < v p+l,q,„ - 2v p,q,n +v p-l,q,n +v p,q+l,n ~ 2v p,q,n + v p,q-l,n >" *>,«?,« ~ zu p,q,n v p,q,n ■ + del 
h 

Assume that /i=& and m,- = = 1, 2. Then simplifying the system to obtain 

h l 

2 

u p,q,n+\ - u p,q,n = m l<" p+l,q,n ~ 2 " p,q,n +u p-\,q,n + "p,q+\,n ~ 2u p,q,n + " p,q-l,n ) " zu p,q,n + zKv p,q,n + zu p,q,n v p,q,n + /* 

2 

l >, 9 ,n+l = m 2< v ;>+!,<?," - 2l >, ? ,« + v p-l,q,n +v p,q+l,n - 2v p,q,n + " M -l,« ' " zKv p,q,n ~ zu p,q,n v p,q,n + del < 

and 

2 

=( 1 - 4m l -2( B+ l))"p, 9 ,„ +m l("^+l, 9 ,« + "p-l,?,n + u p,q-i,n'> + zu p,q,n z + zKv p,q,n + zu p,q,n v p,q,n + & 

2 

v p,q,n+l = (l - 4m 2l v p,q,n + <- v p+l,q,n + u p-\.q,n +u p,q+i,n + " p,q-\,n > m 2 " zKv p,q,n ~ zu p,q,n v p.q.n ■ + del 

This is the alternating direction explicit formula for the Glycolysis model 

Derivation Of Alternating Direction Implicit (ADI)For Glycolysis Model: 

In the ADI approach, the finite difference equations are written in terms of quantities at two x levels. 
However, two different finite difference approximations are used alternately, one to advance the calculations 
from the plane n to a plane n+1, and the second to advance the calculations from (n+7)-plane to the (n+2)- 

d 2 u , d 2 v , 

-^-and — 

or or 



plane by replacing and — ^ by implicit finite difference approximation [5]. we get 



u p,q,n+\ "p,q,n dy dy 

~ ~ 2^ u p+l,q,n+\ ~~ 2u p,q,n+l + u p_iq n+ i) + (u pq + y n - 2lipq n + u pq-\ n ) - u pq n + Kv p q n + 

z h k 

2 

u p.q,n v p.q.n + P> 

v p,q,n+l ~ v p,q,n d^ dj_ 2 

= —( v p+l,q,n+l- 2v p,q,n+l +v p-l,q,n+0 + — ( v p,q+l,n - 2v p,q,n +v p.q-hn*) ~ Kv p,q,n p,q,n v p,q,n + 

z h K 

and 

u p,q,n+2 ~ u p,q,n+\ dy dy 

= _ 2 ("p+lq.n+l - 2u p,q,n+\ + " p-\,q,n+\ ) + ~T ( u p,q+\,n+2 ~ 2u p,q,n+2 + " p,q-\,n+2 ) ~ " p,q,n + Kv p,q,n + 

z h k 

2 

u p.q,n v p.q,n + A 

v p,q,n+2 ~ v p,q,n+l d 2 d 2 2 , , 

= —T <- v p+l,q,n+l ~ 2v p,q,n+l + v p-l,q,n+\ > + ~2 ( v p,q+l,n+2 ~ 2v p,q,n+2 +v p,q-l,n+2>- Kv p,q,n ~ " p,q,n v p,q,n + del 

z h K 
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Simplifying the system will give 

d\z d\z 2 

u p,q,n+l = ~ Y^ u p+l,q,n+l ~ 2u p,q,n+\ + p-\,q,n+0 + ~ 2~'" p,q+\,n ~ 2u p,q,n + u p,q-\,n ) ~ zu p,q,n + zu p,q,n v p,q,n + zKv p,q,n 
h k 

P z + U p,q,n 

d^z d^z 2 

v p,q,n+\ =-^r( v p+l,q,n+l - 2v p,q,n+l + v p-\,q,n+\ ) + (v P,?+l,« ~ 2v p,q,n + v p,?-l,n ) - Kv p,q,n ~ zu p,q,n v p,q,n +zv p,q.n+[ + zdeL 

From the second two equations we will obtain a system of the form 

d\z dp 2 

V9,n+2 = ^(Vl*n+l~ 2 V*n+l + VUn +l) + ^^ 
A it 

A: + » M ,» 

^2^ ^2^ 2 

h k 

From the first two equations where r x = ^y, r 2 = — m x = — 1~ and m 2 = , we get 

/lift k 

"p,q,n+l = r l("p+l,q,n+l - 2u p,q,n+l + " p-l,q,n+l ) + r 2 (" p,q+l,n ~ 2u p,q,n + u p ,q-\n> ~ ™ p ,q,n +zKv p,q,n 
2 

zu p,q,n v p,q,n + P z + U p,q,n^ 
v p,q,n+l = m l( v p+l,q,n+l ~ 2v p,q,n+\ + v p-l,q,n+0 + m 2<- v p,q+l,n ~ 2v p,q,n + v p,q-\,n~> ~ ^ p,q,n ~ 
zu p,q,n v p,q,n + zdel + v p,q.n> 

and this implies that 

2 

u p,q,n+\ = r \ {u p+\,q,n+\- 2u p,q,n+\ + ) + r 2 {u p,q+\,n + " p,q-l,n*> + f 1 ~ 2 1 ~ z >" + zu p,q,n v p,q,n + + zKv p,q,n 

v p,q,n+\ = m \( v p+lq,n+\ ~ 2v p,q,n+l + v p-l,q,n+0 + m 2^p,q+U + V p,q-U ) + (1_ 2m 2 " ^ P ,q,n v p,q,n- + 

We have simplifying these systems of equations yields 

2 

AJ + zKv pqn + zu pqn v pqn 
( l + 2m O v p,q,n+l = m l( v p+l,q,n+l+ v p-l,q,n+0 + m 2( v p,q+l,n + v p,q-l,n^ + ( l - 2m 2) v p,q,n ~ 



*■ p,q,n" p,q,n < 

also the second system of equations, yields 

(1 + 2r 2 )w Pi<?i „ +2 = 1 (« p+l, 9 ,n+l + u p-\,q,n+\ > + < ! " 2 1 )» + r 2 (" p, 9 +l,«+2 + « p, 9 -l,«+2 ) " z " p, ? ,n + 

2 

zu p,q,n v p,q,n + zKv p,q,n + ^ 
(1 + 2m 2 )v pA n+2 = m\ (v p+ l i9i „+l + v p -\, q , n+ \ ) + (1 - 2mj )v p ^ n+ i + m 2 (v p q+1 n+2 + v p q _ l n+2 ) - 
2 

zu p,q,n v p,q,n + zdeL 

The last two systems represent Alternating Direction implicit under the conditions 

u \,q,n+\ = u l,q+l,n+l = 0 
u M,q,n+l = U M ,q+l,n+l = °- 

Also we have 

V l,q,n+l = V l,q+l,n+\ = ^ 
V M ,q,n+\ ~ V M ,q+\,n+\ ~ 

The tri diagonal matrices for the system in the level n advanced to the level n+1, for both u and v can be 
formulated as follows A U=B. 

"(l + 2q) -q 0 0 0.. . 0 

-q (1 + 2IJ) -r x 0 0 0 0 0 0 

0 -r x (l + 2q) -q 0 0 0 0 0 

0 0 -ry (1 + 2))) -/j 

0 0 -I] 

0 0 ... 

0 0 . 



0 0 0 0 0 -jj (1 + 2J]) -ij 

0 0 0 0 0 0 -i) (1 + 2))) 



M 2,g,«+1 
"3,q,n+l 

"5,9+1 



"w-S.g+l 
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& + r 2("2, 9 +l,« + u 2,q-l,n ) + V ~ 2r 2 ~ z >2,q,n + z " 2 ,1," v 2,q,n 



P z + r 2 ("3,q+\,n + U 3,q-Y,n > + ( ! " 



2r 9 



z )"3,<?,« + zu 3,q,n v 3,q,n 



f* + r 2("4,q+\,n + u 4,q-\,n~> + C 1 - 2r 2 ~ z ~> u 4,q,n + 4, 9 ,nv 4 ^ „ 



pz + ^(»M-l, ? +l,/i + "M-1, 3 -1,h) + (! - 2r 2 - Z >"M-1, 9 ,« + zu M -\,q,nv M _ i q n 



(1 + 2mj ) - mj 
- mj (1 + 2mj ) 



0 
0 



(1 + 2mj ) - mj 
- mj (1 + 2mj ) 
0 - mj 
0 0 
0 



0 . . 

0 0 0 
0 0 0 



0 0 -mj a + 2mi) 
0 0 0 -mj 



-mj 

(l + 2»n) 



v 2, 9 ,«+l 
v 3, 9 ,n+l 
v 4,?,n+l 
v 5,?,«+l 



v M-3,^,n+l 
v M-2,q,n+l 
v M-l,q,n+\ 



m 2( v 2,q+l,n + v 2, 9 -l,«) + < J " 2m 2 " ^) v 2, 9 ,n " *» %q,n^2,q,n 

2 

m 2 ( v Xq+\,n + v \q-\,n ) + (1 - 2m 2 - z*>3, 9 ,n - » 3, 9 ,«v 3j9j „ 

2 

m 2 ( v 4,?+l,n + v 4,9-l.n ' + (1 " 2m 2 " zK ^A,q,n ~ zu ^q,n v A,q,n 



m 2( v M-\,q+\,n+ v M-\,q-\,n) + V- 2m 2-? K >M-\,q,n- zu M-\,q, n v M -\,q,n 



And the tridiagonal matrices for the system in level n+l advanced to level n+2 for both u and v are 
given by AU=B. 



(l + 2r 2 ) -r 2 0 0 0.. 

-r 2 (1 + 2*2) -r 2 0 0 0 0 

0 -/j (1 + 2/j) -ij 0 0 0 

0 0 -r 2 (l + 2r 2 ) -r 2 

0 0 -r 2 . 

0 0 

0 0 



0 0 0 -r 2 (l + 2r 2 ) - r 2 
0 0 0 0 -r 2 (l + 2r 2 ) 



"p,2,n+2 
u p,3,n+2 
u pA,n+2 
u p,5,n+2 



u p,M-\n+2 
u p,M-2,n+2 
u p,M~hn+2 
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& + 1 <" p+l,2,n+l + » p-l,2,n+l ) + 2 1 " <■)" p,2,«+l - ZH p,2,n + *» P,2,« v P ,2,/i 

2 

^ + r l(>+l,3,n+l + "H,3,n+l) + ('- 2r r')>,3,»+r z "p,3,n + z " P,l" v p,i,n 

2 

pz + rj (« p+li4iB+ i + « p _i j4 ,„ + i ) + (1-21! - z)n pA „ +1 - zu pA n + zb pA,nv pAn 

2 

^ + l("p+l,5,n+l + " /?-l,5,«+l ) + 0 - 2/^ -z)u p5ill+ i -zu pi5n + zu p,5,nv p ^ n 



P? + r \ <" p+Uf-Vi+1 + » p-l,M-l,n+l ) + (1 - 2/1 - z)h PtM -l,n+l ~ z " pM-\,2,n + z " P^n v pM-l,n+l 



Also for AV=B the tridiagonal is in the form 



(l + 2m 2 ) 


-m 2 


0 


0 


0 


. . . 0 




v p,2,«+2 


-ni2 


(l + 2m 2 ) 


-m 2 


0 


0 


0 0 0 0 




v p,3,«+2 


0 


-m 2 


(l + 2m 2 ) 


-m 2 


0 


0 0 0 0 




v p,4,n+2 


0 


0 


-m 2 


(l + 2m 2 ) 


-m 2 






v p,5,n+2 




0 


0 


-mj 














0 


0 
















0 


0 






















v p,M-3,n+2 


0 


0 


0 


0 


0 


0 -m 2 (l + 2m 2 ) -m 2 




v p,M-2,n+2 


0 


0 


0 


0 


0 


0 0 -m 2 (l + 2m 2 ) 




v p,M-l,n+2 


m l( v p+l,2,n+l 


+ Vl,2,«+1 ) + (1 


-2mj 


-^ v p,2,n+i + zu p,2,n- zu 


2 


P,2,n v p,2,n 



m l ( v p+l,3,«+l + v p-l,3,n+l ) + (1 " 2m l " ^ p,3,n+l + zu p,3,n ~ zu P.3,n v p,3,n 

2 

m l( v >+lA«+l + VW,n+l ) + (1 " 2m l -^) v p,4,«+l + z >,4,» ~ ZM PA<" v p,4,n 



m \ ( v p+l,M-l,n+l + V p-1,M-1,«+1 ) + ( J " 2m l " ^pM-U+l + ^p.M-l.n " PM-l,nv p M _ l n 

2.2. Numerical Stability In Two Dimensional Spaces: 

2.3.1 Numerical stability of ADE: The Von-Neumann method has been used to study the stability analysis of 
Glycolysis model in two dimensions; we can apply this method by substituting the solution in finite difference 

method at time t by y/(t)e m ^ x e mry , where p, y>0 and m = . To apply Von-Neumann on the first 

equation of Glycolysis model 

2 2 

du d u d u 2 

= rfl [ 7T 3 - u + Kv + u v + p, 

dt dx l dy 

2 2 
dv d v d v 2 
— = c/2 [ — T 3 7"] ~ Kv -u v + del, 

dt dx dy 

For the first equation after linearizing it and for some values of p and K neglect the terms pz and 
zKv p q n [3], will be in the form 
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V {t + At)e mPx e m v = (1-4/i - zMt)e m/h e m>y + r x { V {t)e ml3(x+hx) e m7y + 
V {ty^-^) e ^yy + y/ {t)e m ' 3x e mr(y+Ay) + tf/(t)e m/lx e mr(y - Ay) . 
Now dividing both sides of the above equation by lf/(t)e m ^e mn to obtain 

^ (f + A ° = (1 - 4r, - z) + r, (e m/}Ax + e""^ + e m ^ + e^> ) 

= (1 - 4rj - z) + rj (2cos(/?Ax) + 2cos(jAy) 
= (1 - 4/-! - z) + rj(l - 2sin 2 (^Ax/2) + 1 - 2sm 2 (yAy/2j). 
For some values of j8 and y, sill 2 (/?Ax/2) and sin 2 (jAy/2) are unity [4], so 

^>=(l-4,-z) + 2,(-2) 

= (l-4r 1 -z)-4r 1 =(l-8r 1 -z) = <f. 
For stable situation, we need I £ l< 1 , so -1 < (1 - $t\ - z) < 1 , 

2-z 

Case 1: -1 < (1 - 8r, -z) =>r>< 

8 

- z 

Case 2: (1 - 8rj - z) < 1 => q > => neglectthis case because rj > 0. 

8 

For the second (Linearized) equation of Glycolysis model which is in the form 

v p,q,n+l = 0- ~ 4r 2 - Kz)v pqn + r 2 (v p+l q n + v p _ l q n + v p q+l n + v p q _ x n ). 

Let v pqn = (p(t)e m P x e mn and substituting it in the above equation to obtain 

. , , mfk myy ,, . „ . , . mB(x+Ax) myy , , . mB(x-hx) myy 
cp(t + ht)e r e " = (\-\r 2 -Kz)<pit)e 1 e " + r 2 i(pif)e e + 

. . m/& mr( y+Ay) . . m/Jr mr(y-Av). 

Dividing both sides of the above equation by <p{t)e m ^ x e mry to obtain 

(pit + At) mBAx -mBAx myAy -myAy , 

= £ = (1 - 4>2 - Kz) + m 2 [e H + e H + e ' J + e ' J ] 

(pit) 

= (1 - 4m 2 - Kz) + r 2 [2cos(/?Ax) + 2cos(jA)>)] 

= (1 - 4r 2 - Kz) + 2r 2 [1 - 2 sin 2 (/?Ax / 2) + 1 - 2 sin 2 ( jAy / 2)] 

= (1 - 4r 2 sin 2 i/3Ax 1 2) - 4r 2 sin 2 ( jAy /2)-Kz 

= (l-8r 2 - Kz). 

For some values of p and y, we can assume that sin (/?Ax/2)and sin iyAyll) are unity [4], so 

— + = g = (1 - 8^ - Kz) the equation is stable if ||| < 1 which implies that 

(Pit) 

1 1 — 8^2 - Kz l< 1 => — 1 < (1 — 8r 2 — Kz.) < l , which are located in two cases 
Case 1 : - 1 < (l - Sr 2 - Kz) => 8r 2 < 2 - => r 2 < and 

8 

Case 2: (l-8r 2 - Kz) < 1 => 8r 2 > -Kz => r 2 > -Az/8 . 

So the system is stable under the conditions n < - — - ,and r 0 < 2 ~ Kz . 

1 8 2 8 

2.3.2 Numerical stability of ADI: The ADI finite difference form for the first equation of Glycolysis model is 
given in the form 

(1 + 2 ri )u pqn+l = l(Wp+l ;?>n+ l + " p _ [qn+ i) + r 2 iu pq+[n + u pq _ ln ) + 

2 

H-lr 2 -z)u pqn . + zKv pqn + zu pqn v pqn + pz 

2 

By linearization the term zu p q n vanishes, and for the same values of K and p the terms zKv p q n , pz. vanishes. 
Then the equation take the form 



www.ijeijournal.com 



Page I 6 



The Finite Difference Methods and Its Stability for Glycolysis Model in Two Dimensions 



(1 + 2#! )u pqn+i = rj (« p+lqn+1 + u p_l, 9 ,„+l ) + r 2 (« p>9+ i >n + " P;? _l,„ ) + 

In order to study the stability of the above equation, after letting u pqn = y/(t)e m ^ x e mry , where m = 4-1 , we 
obtain 

_ . . . . mBx myy , , . . . m/?(x+Ax:) myy mB(x-Ax) myy . 

/ /.v m7(y+Av) mBx myiy-Ay). , . , . m/& myy 

r 2 iwit)ie e + e 'e ' v ^ y/ ) + (1 - 2rj - z)ys(t)e 'e " . 

Dividing both sides of the above equation by e m P x e m n' to get 

(1 + 21 W* + AO = ^ + A f )( £ m/?Ar + e - m/?Av) )+ ^(0^ _ e "'^ A >') + (1 _ 2q _ zMf) . 
Rearranging, we get 

= ryfit + A0)(2cos(/?Ar) + r 1 (y(f)(2cos(rAy) + (1 + 2r 2 - z)(f(0 

= ^(f + A/) [2(1 - 2sin 2 (/?Ax/2)] + r 2 ^(?)[2(l - 2sin 2 (/?Ay /2)] + (1 - 2q - fc(2> + l))(/(f). 

2 2 

For some values of /? and y, assume that sin (fiAx/2) and sin (j-Ay / 2) are unity [4] , so the equation will take 
the form 

(1 + 2r\)y(t + At) = r lV /(t + Af)(-2) + >y/(f)(-2) + (1 - 2r 2 - z)^(f) 
(1 + 4r 1 )(y(f + Ar) = (1 - 4jj - z)y/(t) 
y/{t + Af) (1 - 4r 2 - z) 



Wit) (1 + 4/!) 



(f(f + A?) 1 - (4r 2 + z) 

So = ^ = £ x . 

V(0 1 + 4#i 

Similarly, for the second equation of the Glycolysis model we assume that v pq n = q>(t)e m ^ x e myy , where 
m = V^i ,the finite difference form of this equation is 

( ! + 2r l >p,q,n+l = r \ (v p+l,q,n+l + v p-l, ? ,n+l ) + r 2 ( V p+1 >9 ,»+1 + v p-l >9 ,n+l ) + 
(l-«z-2r 2 )v p9 _ ln+1 

Then 

. . m/?x mw . ... mB(x+Ax) myy mB(x-Ax) myy. 

(1 + 2r 2 )<p(t + At)e ^ e u = r^<p(t + At)(e ^ e IJ +e ^ v y e ) + 

i m P x myiy+Ay) mBx myiy-Ay). mBx myy 

r 2 i<pit)ie ^ e J + e ^ e ' ) + (1 - Kz - 2^ )<p(t)e ^ e . 

Dividing both sides of the above equation by e m P x e m iy > ( 0 g e ( 

„ . . . . m/fa myy . . . . mB(x+Ax) myy . , . mB(x-Ax) myy . 
(1 + 2r 2 )#>(r + AOe e ' - (p(t + A/)e ^ e + <p(t + At)e ^ e ") = 

i i*\ m P x my(y+Ay) , . m/Jx myy . m& mr(y-Ay). . . m/fe myy 
r 2 {<p(t)e ^ e ' J J -2(p(t)e ^ e " + <p(t)e r e ' J ■" ) + <p(t)e ^ e IJ . 

Which implies that 

(1 + 2r\)cpit + At) = rypit + A0[2cos(yAr)] + r 2 <p(0[2cos(yAy)] + (1 - Kz - 2r^)<p(t) 

(1 + 2r x )<p(t + AO = ryp(t + AO [2(1 - 2sin 2 (/?Ax/2)] + r 2 (2(l - 2 sin 2 (yAyl 2)<p(t) + (l-Kz- 2r )p(f). 

For some values of p and y, we have sin (/?Ar/2)and sin (jAy/2) are unity, so we have 

(1 + 2^)^ + AO = r x q>{t + Af)(-2) + r 2 ^(0 + (l-Kz- 2r )(p(t). 

(1 + 4q )<p(t + At) = (1 - Kz - 4r 2 )(p(t) 

y/it + At) _ (1 - (4r 2 + Kz) _ 
Wit) (1 + 4^) 2 ' 

Where ^ and ^ stand for the I-plane and Il-plane respectively, each of the above terms is conditionally 

stable. However the combined two-levels has the form 

(1 - (4r 2 + Kz)) (1 - (4r 2 + Kz)) 

s =££ =r- — ir — — — — i 

\/)/ 12 : !l 



(1 + 4^) (1 + 41!) 
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Thus the above scheme is unconditionally stable, each individual equation is conditionally stable by itself, and 
combined two-level is completely stable. 

III. Application (Numerical Example) 

We solved the following example numerically to illustrate the efficiency of the presented methods, 
suppose we have the system 

— = di Au — u + kv + u 2 v + p, 

dt 1 r 

— = d 7 Av — kv — u 2 v + S 

dt L 

We the initial conditions 

U(x, 0)=Us+ 0.01 sm(nx/L) f 

or 0<x<Z 

V (x, 0) = Vs - 0. 12 sin(nx/L) for 0 <x < L 

U(0, t) = Us, U(L, t) = Us and V (0, t) = Vs, V (L, t) = Vs 

We will take 

d 1= d 2 =0.01 , p=0.09 ,5=-0.004 ,U S =0 , V s =l 
Then the results in more details are shown in following table and figure: 



1.0003 


1.0000 




1.0017 


1.0013 


1.0045 r 


1.0028 


1.0024 




1.0036 


1.0032 


1.004 


1.0040 


1.0038 


1.0035 


1.0042 


1.0041 


1.0042 


1.0043 


1.003 


1.0041 


1.0044 


1.0040 


1.0044 


1.0025 


1.0038 


1.0044 




1.0038 


1.0044 


1.002 


1.0038 


1.0044 




1.0040 


1.0044 


1.0015 


1.0041 


1.0044 


1.001 


1.0042 


1.0043 


1.0042 


1.0041 


1.0005 


1.0040 


1.0038 


1.0036 


1.0032 


1 

D 


1.0028 


1.0024 


1.0017 


1.0013 




1.0003 


1.0000 





u(K.t) 





20 



IV. 



Conclusion 



We saw that alternating direction implicit is more accurate than alternating direction explicit method 
for solving Glycolysis model and we found that ADE method is stable under condition r r < — -, and < — — , 

8 8 

while ADI is unconditionally stable. 
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